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ABSTRACT 

Based on new calculations we reconfirm the low and high density limits on the for- 
bidden fine structure line ratio [Oil] 7(3729)//(3726): limjv ^o = 1-5 and limjv,,-,^ 
= 0.35. Employing [O n] collision strengths calculated using the Breit-Pauli R-matrix 
method we rule out any significant deviation due to relativistic effects from these 
canonical values. The present results are in substantial agreement with older calcula- 
tions by Pradhan (1976) and validate the extensive observational analysis of gaseous 
nebulae by Copetti and Writzel (2002) and Wang etal. (2004) that reach the same 
conclusions. The present theoretical results and the recent observational analyses dif- 
fer significantly from the calculations by MacLaughlin and Bell (1998) and Keenan 
etal. (1999). The new maxwellian averaged effective collision strengths arc presented 
for the 10 transitions among the first 5 levels to enable computations of [O n] line 
ratios. 
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1 INTRODUCTION 

The most prominent density diagnostics in astrophysics is 
due to the forbidden fine structure lines [O n] AA 3729, 3726 
and [S n] AA 6716, 6731. Their utility stems from several fac- 
tors such as: (A) they respectively lie at the blue and the red 
ends of the optical spectrum, (B) their atomic structure and 
hence the density dependence is essentially the same, and 
(C) they are quite strong in the spectra of most H II regions 
owing to the relatively large abundances of oxygen and sul- 
phur. Seaton and Osterbrock (1957) have described the basic 
physics of these forbidden transitions. High-accuracy calcu- 
lations using the then newly developed computer programs 
based on the close-coupling method IMPACT (Eissner and 
Seaton 1972, 1974, Crees etal. 1978) were later carried out 
by Pradhan (1976; hereafter P76) for the collision strengths, 
and by Eissner and Zeippen (1981) and Zeippen (1982) for 
the transition probablities. These atomic parameters subse- 
quently enabled a consistent derivation of electron densities 
from observations of [O n] and [S n] lines in a wide variety 
of H II regions (e. g. Kingsburgh & English 1992, Aller & 
Hyung 1995, Aller etal. 1996, Keyes etal. 1990). 

More recently, McLaughlin and Bell (1998; hereafter 
MB98) repeated the [O n] calculations of collision strengths 
using the R-matrix method (Burke etal. 1971; Berrington 
etal. 1995), also based on the close coupling approximation 
and widely employed for a large number of atomic calcula- 
tions {The Opacity Project Team 1995, Hummer etal. 1993). 
They included a much larger target wavefunction expansion 



than P76, and relativistic effects not considered in the ear- 
lier calculations. Their electron impact collision strengths 
and rate coefficients were markediy different for the rele- 
vant transitions 4 S3/ 2 — >2 Dg / / 2 3 , 2 than P76, which leads to 
the theoretical density diagnostic line ratio 7(3729)//(3726) 
to be up to 30% higher and « 2.0 in the low-density limit 
limjv c — > 0. Keenan etal. (1999; hereafter K99) recomputed 
the [O n] line ratios to analyze several planetary nebulae 
using these MB98 results. 

However, other extensive observational studies (e. g. 
Copetti and Writzel 2002, Wang et al. 2004) have noted the 
discrepancy between electron densities derived from [O n] 
and other density indicators, notably the [S n] AA 6716, 6731. 
In particular the recent analysis of a sample of over a hun- 
dred nebulae by Wang et al. (2004) shows that the colli- 
sion strengths of MB98 are not supported by observations, 
and that the earlier results of P76 are to be preferred. But 
these observational studies leave open the question of what 
precisely are the collision strengths. Given that P76 used 
a small basis set to describe the O II target, considered no 
relativistic effects, and could not fully resolve the resonance 
structures in the collision strengths owing to computational 
constraints, it seems puzzling that the new MB98 results 
which do account for all of these factors appear to be inccu- 
rate. 

To address this important issue and to resolve the out- 
standing discrepancy, we recently undertook new calcultions 
for [O n] using the same Breit-Pauli R-matrix method as em- 
ployed by MB98 and including relativistic effects. Whiie the 
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details of the atomic calculations and comparison of colli- 
sion strengths with different basis sets of target wavefunc- 
tions will be presented elsewhere (Montenegro etal. 2005), 
we present the final results of astrophysical interest in this 
Letter. 



2 THEORY AND COMPUTATIONS 

Forbidden lines are often sensitive to ambient electron den- 
sity; as the Einstein spontaneous decay rates of the upper 
levels are small, they may be collisionally excited to other 
levels by electron impact before radiative decay (Osterbrock 
1989, Dopita and Sutherland 2003). This is likely to happen 
when there is a pair of lines originating from closed spaced 
metastable energy levels, especially in ions of the 2p 3 and the 
3p 3 outer electronic configurations as exemplified by O II and 
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etal. 2005). In the present calculations we employ a 



S 11. The first five levels are: 4 S 3/2 , 2 Dg/ 2 3 , 2 , 
pair of [On] transitions of interest are 2 Dg / / 2 3/2 
at A 3729 and A 3726 respectively. 

The basic physics of the limiting values of the line ratio 
7(3729)/7(3726) is quite simple. At low electron densities 
every excitation to the two metastable levels 2 D° 5/ / 2 3 / 2 ) is 
followed by a spontaneous decay back to the ground level 
4 S 3 / 2 since the collisional mixing rate among the two excited 
levels is negligible. In that case the line ratio in principle 
must be equal to the ratio of the excitation rate coefficients 



lim 



1(3729) g( 4 SS /2 - 2 Dg /a ) 



_2 ]-)o \ 

/2 U 3/2> 



iVc-o 1(3726) q( 4 S° /2 
where the excitation rate coefficient qtj is 



.63 x 10" 6 cm 3 /s ■ exp{-Eij/kT e 



9 



Ty(T.; 



(1) 



(2) 



where gi is the statistical weight of the initial level and the 
quantity Ty is the Maxwellian averaged collision strength: 



Q ij (E)exp(-E/kT e ) d(E/kT e ). 



(3) 



If relativistic effects are negligible then the collision 
strengths may be calculated in LS coupling, and an alge- 
braic transformation may be employed to obtain the fine 
structure collision strengths. This was the procedure em- 
ployed in P76. The ratio of fine structure LSJ to LS col- 
lision collision is especially simple when the lower level has 
either L or S = 0, such as for O 11 and S 11, i. e. 



Q(SLJ - S'L'J') 
Q(SL - S'U) 



2J' + 1 



(2S' + 1)(2L' + 1) 



(4) 



If the excited levels are so closely spaced that the excita- 
tion rates have virtually the same temperature dependence, 
the line ratio would then be equal to the ratio of the sta- 
tistical weights (2J'+1) of the upper levels 2 Dgy 2 >3 / 2 , i.e. 
6/4. If, however, relativistic mixing is significant then the 
line ratio will depart from the LS coupling value. That was 
the contention of MB98 and K99. 

Therefore we carry out the present calculations 
including relativistic effects and with a suitably large 
target wavefunction expansion. The Breit-Pauli R-matrix 
(BPRM) calculations with different trial basis sets 
are described in detail in another paper (Montenegro 



16-level target: ls 2 2s 2 [2p J ( 4 S^ /2 , J D° /2i3/2/ ,^P° /2il/2 ); 

2s2p^ ( P5/2, 3/2, 1/2)^-^5/2, 3/2)5 

2p 2 3s( 4 P lA3A5/2 , 2 P 1/2 ,' 3/2 ), 2s2p 4 ( 4 S 1/2 )]. This close 
coupling expansion is more than sufficient to obtain 
accurate collision strengths for the first 5 levels. All 
resonance structures up to the highest target threshold 
energy £(2s2p 4 ( 4 Si /2 )) = 1.7829 Ry are resolved. The last 
threshold lies sufficiently high to ensure that resonance and 
coupling effects in the collision strengths are fully accounted 
for in the excitation of the first 5 levels considered in the 
collisional-radiative model to compute the line ratios. 
Collision strengths at energies 1 Ry higher than the highest 
of the first 5 levels {E( 2 P° /2 ) = 0.369 Ry) contribute 
negligibly to the rate coefficients; at T = 20,000 K the 
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accordingly for E > lRy in Eq. (3). 

In the calculation of [O 11] line ratios we employ the 
transition probabilities from Zeippen (1982). Eissner and 
Zeippen (1981) computed the A-values for [O 11] transitions 
taking full account of the magnetic dipole Ml operator, and 
they showed that in the high density limit the line ratio 
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3 RESULTS AND DISCUSSION 

Fig. shows the fine structure BPRM collision strengths 
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These figures appear to be the first clear presentation of the 
resonances in these collision strengths. P76 did not present 
detailed resonance structures except in the near-threshold 
region of n( 2 Dg /2 - 2 D 3/2 ). MB98 plotted these on an 
energy scale up to 5Ry, well above the resonance region 
up to ~ 2 Ry, but which does not exhibit the resonances 
in detail to enable comparison. An interesting feature clear 
from Fig. Q is that the resonances do not play a large 
role in ^( 4 S 3/2 - 2 D° /2 ) and n( 4 S^ /2 - 2 D° /2 ) and hence 
the rate coefficients for these transitions. Although they 
are significant in Q.( 2 T)°^ 2 — 2 D 3j/2 ), collisional mixing via 
this transition is not important in the low density limit, 
which therefore depends only on the ratio of the transitions 
from the ground state 4 S 3/ / 2 to the 2 Dg / / 2 3 , 2 levels. We 
find that this ratio is constant at 6/4 throughout the 
energy range under consideration, resonant or non-resonant 
values. Therefore we do not find any signifant evidence 
of relativistic effects, which would manifest itself in a 
departure from this ratio. Remarkably the present total 
sum ^2 j = 5/ 2 3/2 ^( 4 ^3/2 — 2 D j) ~ 1.42, compared to the 
LS coupling P76 value of 1.31, and an even earlier value of 
1.36 obtained by Saraph, Seaton, and Shemming (1969). 

Table 1 give the Maxwellian averaged collision strengths 
T(T) for the 5-level [On] model. At 10,000 K we obtain 
T( 4 S 3/ / 2 — 2 D 3/ / 2 ) = 0.585, in good agreement with the earlier 
P76 value of 0.534, about 9 % lower, but considerably higher 
than the M98 value of 0.422 (quoted in K99) which is 28% 
lower than the new value. More importantly our results dis- 
agree with MB98 for the ratio discussed above. It is this ratio 
that is responsible for the K99 line ratio 7(3729)/7(3726) to 
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Figure 1. Collision strengths for the fine structure transitions associated with the [Oil] line ratio 3729A/3726A. 
Note that fi( 4 S°y 2 — 2 D°^ 2 )/Q( 4 S°y 2 — 2 D°^ 2 ) = 1.5 througout. There is significant resonance enhancement in the 

collisional mixing transition 2 D5/ 2 — 2 ^3/2> mse * shows the near-threhold resonances on an expanded scale. 



Table 1. Effective Maxwcllian averaged collision strengths 
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0.151 


0.152 
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0.156 


0.158 


0.383 


0.392 


0.402 


0.414 


0.428 


0.441 


0.376 


0.386 


0.397 


0.409 


0.423 


0.437 


0.277 


0.284 


0.291 


0.300 


0.310 


0.321 



be about 30% higher (~ 2.0) than the expected low-density 
limit of 1.5, as shown in Fig. 2. 

Comparing the present relativistic BPRM results for 
effective collision strengths with the LS coupling results of 
P76 we find good agreeement, mostly within a few percent, 
with the notable exception of T( 2 D5/ 2 — 2 D^/ 2 )- Owing to 
the more extensive delineation of resonance structures in the 
present calculations (Fig. 1), the T value is much higher. For 
example, at T = 10,000 K the P76 value is 1.168 compared 
to the present value of 1.426 in Table 1. On the other hand 
the present T( 2 P3 /2 - 2 P° /2 ) = 0.291 agrees well with the 
P76 value 0.287 at T = 10,000 K. 



Fig. 2 shows that the present collision strengths yield 
the line ratio Ri = J(3729)/I(3726), which approaches 
the low- and high-density limits exactly. The difference is 
not discernible when we use the P76 values. On the other 
hand the difference with MB98 is quite pronounced and 
approaches ~30% in the low-density limit. The tempera- 
ture variation between T = 10,000 K (solid line) and T = 
20,000 K (dashed line) is also small, demonstrating the effi- 
cacy of this ratio as an excellent density diagnostic. 
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Figure 2. [O II] line ratio I (3729)// (3726) vs electron density N e : 
present results — solid line, Pradhan (1976) -.- dot-dashed line 
(nearly indistinguishable from the solid line), McLaughlin and 
Bell (1998) .... dotted line, at T = 10,000 K. The dashed line is 
the line ratio at T = 20,000 K. 

4 CONCLUSION 

We have carried out new relativistic Breit-Pauli R-matrix 
calculations for the [O n] transitions responsible for the im- 
portant density diagnostic line ratio 3729A/3726A. We find 
no evidence of any significant departure from the earlier LS 
coupling results of Pradhan (1976). The line ratios derived 
from the present results also agree with the canonical limits 
expected on physical grounds. The new results are in con- 
siderable disagreement with the calculations of McLaughlin 
and Bell (1989) and the line ratios of Keenan etal. (1999). 
We also reconfirm the observational analyses of Copetti and 
Writzel (2002) and Wang etal. (2004). 
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